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ABSTRACT 


A reinvestigation was conducted of an anomaly in the experimental 
angular correlation for the gamma-gamma cascade of Ni®9 which was 
reported by Verser and Conway in 1960 and 1961, respectively. The measure- 
ment of this correlation was being used to test angular correlation 
equipment built at the US Naval Postgraduate School. Two new experimental 
determinations of the correlation function with the same equipment as used 
previously indicated good agreement with the theoretical function. 
Verser's and Conway's experimentally-determined angular correlation 
functions were found to deviate less significantly from the theoretical 
function than they had reported since their standard deviations had failed 
to take into account the variations of the characteristics of the electronic 
equipment due solely to electronic causes. It is concluded that the anomaly 
reported by previous experimenters was statistical in nature. 


Programs for the CDC 1604 computer have been developed in FORTRAN 
language for complete data reduction accepting as input the total counts 
in each channel, the total coincident counts, and other experimental 
Parameters. The program provides for the calculation of the coefficients 
for a cosine series and a Legendre polynomial series expansion of the 
least squares fit to the experimental points. Other programs were 
developed to calculate the corrections for non-point detectors to be 
applied to the theoretical coefficients, and to calculate coincidence 
resolving times from data obtained during a large number of runs. 


Major improvements were made in the equipment to be used in future 
angular correlation experiments, including the construction and testing 
of both a slow and a fast coincidence unit. 


The author wishes to express his appreciation to Associate Professor 
H. E. Handler who has guided the development of angular correlation 
measuring equipment at the US Naval Postgraduate School for the last 
three academic years. His recommendations, particularly in the area 
of data analysis, and his perusal of the draft of this paper resulted in 
a still more thorough reinvestigation of the results of the angular 
correlation studies reported by previous experimenters at this School. 
It is a pleasure to ackowledge the guidance and aid received from Lt Cmdr 
P, A. Phelps (USN) in comprehending the theoretical aspects of the basic 
experiment, in becoming familiar with the electronic equipment, and 
in writing and debugging the computer programs. Appreciation is expressed 
to Mr. M. A. Brillhart for his construction of the coincidence units, and 
for his assistance in developing the mixer circuitry. 


ie 





TABLE OF CONTENTS 





Section Tike Page 
i ie introduction 
A. Purpose of Angular Correlation Experiments 1 
B. Methods of Conducting the Experiment 2 
C. Interpretation of the Experimental Data 5 
D. Previous Work in this Field at USNPS 4 
E. Scope of the Investigation 6 
TT. Reinvestigation of the Earlier Reported Anomaly 
A. Equipment Used 8 
B. Coincidence Resolving Time Tests 8 
C. Daily Checkout Procedure and Instrument log Pt 
D. Scaler Dead Time ll 
E. Calculation of the Standard Deviation of W694) 1 
F. Preparation and Check of New Computer Programs 16 
G. Reanalysis of Earlier Data 18 
H. Analysis of New Data 75S) 
I. Comparison of Experimental Results with Theory Ze 
ET Development of New Instrumentation 
A. New Slow Coincidence Unit ay} 
B. Requirement for a Fast Coincidence Unit zy 
1. General 27 
2. Relation Between r{w..@,)| and TC.) on 
3. BELect on o-(6 ) of Decreasing the 33 
Coincidence Resolving Time 
4. Relative Effect of Coincidence Resolving 59 
Time Drift 
C. Operation of the Fast Coincidence Unit 34 
D. Testing and Alignment of the Fast 36 
Coincidence Unit 
E. Equipment Status 42 


Appendix FORTRAN Programs 





A. Program for the Least Squares Analysis 44 
B. Program for the Solid Angle Correction 60 
C. Program for the Coincidence Resolving Time 66 


ak at 








=_— ao = = 


— < 
ae r 





— 
=) i 


=.= 





—=_ 
a 





OO 





= 








D. Results of Analysis by the Least Squares 69 
Program 


E. Glossary of Terms for the Least Squares Analysis 7/7 


Bibliography 19 








LV 

















Figure 


i 


oO ON DO F&F Wb 


Re RP FF 
WwW Nh F&F © 


14. 
, 
16. 


LIST OF FIGURES 


Title 
Decay Scheme of cae 
Results of Verser and Conway 
Experimental Setup of Verser and Conway 
Angular Correlation Table 
Instrument Log 
Reanalysis of Data of Verser and Conway 
Analysis of New Data 
Composite Data Analysis 
Block Diagram of New Slow Coincidence Unit 
Circuit Diagram of New Slow Coincidence Unit 
Block Diagram of Experimental Setup 
Block Diagram of Fast Coincidence Unit 


Fast Coincidence Pulse Amplifying and Discriminating 
Circuitry 


Circuit Diagram of Mixer Circuit 
Fast Coincidence Unit Test Scheme 


Electronic System 


Page 


10 
12 
19 
24 
26 
28 
29 , 30 
30 
35 
37 


38 
40 


43 





LIST OF TABLES 





Table Title Page 

Results of Reanalysis of Earlier Data ZG 

2) Effect on W, (9) of Varying Input Parameters 20 

oe Effect on Solid Angle Corrections of Varying ap) 
Input Parameters 

4, Results of Analysis of New Data Ze 

or Results of Composite Data Analysis Zo 




















vi 














INTRODUCTION 
Purpose of Angular Correiation Experiments 


In gamma-gamma angular correlation experiments, one observes 
the probability of emission of the second quantum of a cascade as 
a function of angle relative to the first quantum. In general, 
this probability is an anisotropic function of the angle. Such 
experiments, in conjunction with other nuclear studies, can lead 
to identification of the multipole orders of the emissions and 
the total angular momenta of the original, intermediate, and final 
states of the nuclide under investigation. j1,2,3 | The orienta- 
tion of the total angular momentum vector of the intermediate level 
is correlated with the direction of emission of each quantum. The 
anisotropy of the probability distribution pattern will be attenuated 
if the orientation of the total angular momentum vector of the inter- 
mediate level changes between the times of these two emissions. Thus, 
one must study a nuclide for which the intermediate energy level has 
a mean life which is short relative to the period of precession of 
the angular momentum vector of that level. Only if the total angu- 
lar momentum vector of the intermediate level does not appreciably 
change its orientation between the emissions of the two quanta can 
the probability distribution as a function of angle between the two 
emissions be interpreted meaningfully. 

Generally, no parity assignments can be made for the nuclear 
levels involved based solely on a directional angular correlation 
study such as described above, since the angular distribution of 
the emissions is the same for electric and magnetic multipole radi- 
ations of the same order. By means of a similar polarization- 
direction correlation experiment, wherein one measures the orientation 
of the polarization vector relative to the plane of the two 
quanta, one can determine whether the emissions are electric or mag- 
netic multipoles. [4,5] This information, coupled with known mul- 
tipole orders obtained from an angular correlation study, enables 


the assignment of relative parities to the three nuclear energy 
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levels. Absolute parity assignments require additional data. 
In the case of wie? an even-even nuclide, an even parity 
and zero total angular momentum are predicted for the ground 
state on the basis of the nuclear shell model, and are corrobo- 
rated by other experimental techniques unrelated to this paper. 
Studies of the decay scheme have shown that the spins of the 
first and second excited levels are 2 and 4 respectively. Cor- 
relation studies are consistent with this and demonstrate that 
these radiations are pure electric quadrupoles. Thus the com- 


plete decay scheme for vee is 
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Since this decay scheme is well understood, this gamma-gamma 
cascade of io? is frequently used to determine whether angular 
correlation equipment is functioning properly. | 6| It is used 


in this manner in this paper. 
Methods of Conducting the Experiment 


Since nuclei are generally randomly oriented, one cannot 
measure the radiation pattern of a given transition. One way of 
minimizing the random orientation of the nuclei is to apply a 
strong magnetic field or a strong electric field gradient at very 
low temperatures; radiation patterns have been measured using this 
technique. In the event of a cascade, the effects of this random 
orientation can be partially circumvented by establishing some re- 
ference direction with respect to which only certain of the second 


of the two radiations are measured. In the experiment discussed in 





this paper, nuclei have not been aligned; the reference direction 
for observation is the direction of emission of the first quantum 
in the two-gamma cascade of re Since the lifetime of the inter~- 
mediate nuclear level for this nuclide is short in comparison with 
the period of precession of its angular momentum vector, the ori- 
entation of the angular momentum vector does not vary significantly 
between the times of emission of the two quanta. Consequently, an 
unattenuated angular correlation between the directions of emission 


of the two successive radiations is observed. 


Interpretation of the Experimental Data 

By measuring the total,coincidence counting rate at several 
angles 9 between successive quanta, and applying corrections for 
accidental coincidences and for the attenuation of the correlation 
function due to the use of detectors of finite size, 7,8] one can 
determine the relative likelihood of emission of the second quantum 
at these angles relative to the direction of emission of the first. 


The experimental estimate of the correlation function at a specific 


angle Oo. is calculated according to the formula 
C 
W. (8) =|—2— (1) 
ex’ i N 
| ae 


where C. is the genuine coincidence counting rate at the given angle, 
and Ny and Ny are the individual channel counting rates at that 
angle corrected for the system dead time and background. The least 
Squares fit to these data is the experimental correlation function 
W, 69) - 17] W, 69) is expressed in terms of a linear combination 

of even powers of cos 8 or in terms of a linear combination of even 
order Legendre polynomials. In the case of the gamma-gamma cascade 
ron we the theoretical correlation function, W469) » corrected for 
the solid angles -subtended by the detectors used in this experiment 


is given by 


W 69) = 1.00000 + 0.09960 P,(cos 0)+ 0.00838 P, (cos 0) . (2) 














Experimental results can be compared with theory by comparing the 
coefficients calculated in the least squares fit with the above 
coefficients. In addition, one can compare the experimentally- 
determined anisotropy with the theoretical value. The anisotropy 


is defined as follows: 


O 
W, , (180°) 


= 7 (900) ee (3) 
W, 690°) 


R 


The theoretical value of R, corrected for the geometry of this 
experiment is + 0.16221. A close agreement between the experi- 
mental and theoretical coefficients and between the experimental 
and theoretical anisotropies provides a check of proper equipment 


operation. 


Previous Work in This Field at U.S. Naval Postgraduate School 

The basic equipment for this experiment was assembled in 
academic year 1959-60 by Verser be who conducted an angular 
correlation experiment in the 90-180° quadrant and reported agree- 
ment with the theoretical anisotropy but reported marked deviation 
from theory in the shape of the experimental curve. The equipment 
was modified slightly in academic year 1960-61 by Conway [10 } who 
performed the same experiment in the same quadrant and also reported 
a similar deviation from theory in the shape of the experimental 
curve. Fig (2) shows the results of these investigations. Verser 
reported a dip in his experimental curve which deviated by more than 
two standard deviations from the theoretical curve in the 120-135° 
region. Conway reported a similar dip in the vicinity of 135). and 
a less pronounced dip in the 120 region. 

Conway wrote two programs in CDC 1604 computer assembly routine 
language, one for the least squares fit of a linear combination of 
Legendre polynomials of even order to the experimental data, and a 
second to calculate the solid angle corrections necessitated by the 
finite size of the detectors. The least squares program required a 
considerable amount of hand-calculated data reduction before it could 


be used. 








OBT SO Osi 





O°T 
(*@) S1IO}JIeI0q UseMJeq sTsuy __ 
“ e a oe 
ST 
ABMUO) aoe! ° a 
I ™ . i 
(0) y . 
: | 
i oe < | 
IBSADA a ¢ / | 
("e) M | 
| 4 “\s0°T 
‘ i 
/ / | 
p | 
a | ur) 
ad : 
Ve 
4 . 
/ ee COL i ay 
(9) M EE 
\f ' 
A Y AB®Mu0) (06) M4 Lepomis 
/ (0) iM 


VY (06) Mm 


A@MuoD pue AasIdsA Jo S A[Nssy 


Zi T# PANST A 


al 
(9)¢ M 
oa 























Possible mechanisms underlying the reported anomalous dip 
in the experimental curve as suggested by Verser and Conway in~ 
clude the following: 

1. A change in the conditimof the source, possibly 
crystallization resulting in preferential align- 
ment of the nuclei within the crystal due to in- 
ternal electric fields. 

2. Fluctuations of ambient temperature which resulted 
in changes in the operating characteristics of the 
equipment. 

3. Uncompensated variations in line voltage. 

4. Perturbing effects of external magnetic fields. 

5. Improper orientation of the source relative to the 
detectors. 

It was recognized that the likelihood of having any noticeable 
effects due to Items "1" or "4" above at room temperature was ex- 
tremely small since the Boltzmann factor for the relative popula- 
tions of the energy levels is nearly unity except under the influ- 
ence of extremely strong fields. 

In addition, the possibility was recognized of there being in 
fact no real anomaly, but that the discrepancy was simply an un- 
likely statistical occurrence or that some systematic error existed 
in the experimental technique, in the processing of data by hand, or 
in one of the previous computer programs. The likelihood of this 
deviation from the theoretical being simply a statistical fluke be- 
came greater when it was recognized that the variation in the re- 
solving time of the coincidence unit was not taken into account in 
computation of the standard deviation of the correlation function 


evaluated at the angles at which the experiments were conducted. 


Scope of the Investigation 

Many experimenters have reported agreement with the theoretical 
gamma-gamma angular correlation function for wie? [12.12]. yet two 
experimenters at this School independently reported a significant 











deviation from theory in the 105-135° region. Thus it was clear 
that before one could use this equipment for the experiments for 
which it was assembled; namely, experiments with short-lived nu= 
clides, further investigation was necessary either to arrive at 
a mechanism for the reported anomaly, or to demonstrate that the 
anomaly, under more closely controlled experimental conditions, 
is not observed. In the investigation of this anomaly, certain 
improvements were made in the sepeninenead equipment; new pro- 
grams in FORTRAN language were developed to compute the solid 
angle corrections to the angular correlation function and cal- 
culate the least squares fit; and experiments were conducted 
under closely controlled conditions in both the 90-180° and the 
180-270" quadrants to see whether a similar dip would be observed. 
Two additional coincidence units were constructed. One of 
these is a slow coincidence unit which, unlike that used by Verser 
and Conway, has a resolving time which is independent of counting 
rate. This characteristic will be required in conjunction with 
future experiments with short-lived nuclides. The second coinci- 
dence unit has a resolving time an order of magnitude shorter than 
that of the slow units. The use of the fast unit in conjunction 
with a slow unit enables one to achieve better precision by reducing 


the accidental coincidence counting rate. 








103 
REINVESTIGATION OF THE EARLIER REPORTED ANOMALY 


Equipment Used 

The experimental setup employed by Verser and Conway is des- 
cribed in detail elsewhere [9,10]. Preune (3) indicates the general 
test setup and Figure (4) shows the angular correlation table with 
the source, crystals, photomultiplier tubes, and preamplifiers. The 
same slow coincidence unit as used by Verser and Conway was used in 
the angular correlation experiments in the 90-180° and in the 180- 
270° quadrants. The only significant change in the setup for these 
experiments consisted of adding a voltage regulator to minimize the 
effects of uncompensated line voltage fluctuations. Regulated vol- 
tage was then supplied to those parts of the electronic apparatus 
most likely to distort the experimental measurements when under the 
influence of line voltage fluctuations. These included the coinci- 
dence unit, both linear amplifiers, and the high voltage power supply 


for the photomultiplier tubes. 


Coincidence Resolving Time Tests 

The accidental coincidence counting rate is calculated from the 
usual formula 

C, = 27 N,N), (4) 

where 7 is the coincidence resolving time. Before and during the 
angular correlation experiments, scores of tests were conducted to 
determine the stability and the magnitude of 7 . The unit exhibited 
random fluctuations in resolving time as much as three times as large 
as those that can be explained in terms of nuclear statistics. The 
fractional variation in T , a was observed to be approximately 
twice that reported by Verser and Conway. The measured mean is 
0.2470 microseconds. The variance of a single measurement of T is 
0.0044 microseconds; this includes the variation due to non-nuclear 
fluctuations and was used in calculating the variance of the acci- 


dental coincidence counting rates in this work. Since the resolving- 





Figure 3. Experimental Setup of Verser and Conway 
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Preamp-2 Shaper-2 Amp-2 
Channel 2 PHS~2 : Scaler~-2 | 


Xtal's = Harshaw Scintillator Crystals, Type &$8/2.(2" diam x2" thick) 


PM's = RCA 6342 Photomultipler Tubes 
Power for PM tubes from Hamner H.V. Power Supply Model N-~401 








Amp's = Hamner N302 Amplifiers - each includes a differential pulse 
height analyzer which supplies a standard pulse out (''PHS" out) 
(A linearly amplified input is also supplied as an output but 
was not used in this experimental setup). 





Coinc = Coincidence Anzulyzer, Atomic Model 502A. 


Scalers = Baird-Atomic, Model 1283 
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time characteristics were investigated and documented far more com~- 
pletely than in previous years, the observed instability was felt 
to be similar to that which would have been found by the earlier 
experimenters had they investigated this problem more thoroughly. 
An inclusion of the effect of this fluctuation in the analysis of 


their data would decrease their reported precision. 


Daily Checkout Procedure and Instrument Log 

A daily checkout procedure was devised to ensure that exactly 
the same checks of the operating characteristics of the electronic 
equipment were performed each day in the same sequence and with the 
equipment in the same operating configuration at each stage of the 
checkout. 

» An instrument log was developed in which the experimental 
equipment layout was recorded, and monitoring points for key wave 
shapes were indicated. Photographs of these key wave shapes are 
included, as well as the record of the oscilloscope settings used 
when the photographs were taken; a sample page from this log is 
shown in Figure 5 . This log has proven useful in trouble-shooting 
the extensive circuitry and for monitoring changes in operating 
characteristics, such as pulse rise times, which could degrade system 


performance. 


Scaler Dead Time 

As a part of the complete checkout and alignment of all elec- 
tronic gear prior to conducting the angular correlation experiments, 
two-source experiments were run to determine the system resolving 
time (primarily scaler dead time). Several of these measurements all 
led to the same surprising result, namely, a negative dead time for 
the system. Further investigation of this problem involved supplying 
pairs of pulses from a double pulse generator at a constant repetition 
rate and varying the time between the pulses of each pair. Expected 
behavior was that with the 1000 per second repetition rate for the pair 
of pulses, 2000 counts per second would be registered on the scalers as 


long as the separation between the pulses of each pair exceeded the 
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resolving time of the system. When the separation between the two 
pulses was less than the resolving time of the system, 1000 counts 
per second would be expected. What was observed was the recording 
of 2000 counts per second when the interval between the pulses of 
each pair was greater than about 5-10 microseconds, then as the in- 
terval was decreased, the count rate recoracampy tie seavers in; 
creased, on occasion by as much as 25%; i.e., to about 2500 counts 
per second. 

Since the experimentally-determined correlation function is 
not sensitive to the introduction of these spurious counts because 
the error introduced into the single channel counting rates is es- 
sentially constant for all angles and the effect on the coincidence 
counting rate is negligible, this anomalous behavior had no signifi- 
cant influence on the experimental results reported herein. New 


scalers will be employed for future experiments. 


Calculation of the Standard Deviation of W694) 

The method of calculating the standard deviation of the experi- 
mental estimate of the correlation function has been the subject of 
extensive consideration. The experimental estimate of the correlation 


function at a specific angle oO. is calculated according to the formula 


e 





| Co 
Ww (0.) = f= (1) 
et 
where C. is the genuine coincidence counting rate at the given angle, 
and Ny and Ny are the individual channel counting rates at that 


angle corrected for the system dead time and background. The genuine 


coincidence counting rate is calculated from 


c,=C,-C 5 

G ay As (>) 

where Co = the genuine coincidence counting rate, Cy. = the total coinci- 

dence counting rate, and Cy = the accidental coincidence counting rate. 
13 











C.. and the single channel counts were measured directly by scalers 
and the accidental coincidence rate was calculated from 

C= 2TN, Ny. (6) 
It can be seen that the standard deviation of Wax Op? depends 
almost entirely upon that of Co: The standard deviation of C is 
obtained by taking the square root of the sum of the variances of 


the two presumably independent quantities, on and C, 


7 (C.) = r2¢6,) BO Ge. (7) 


Verser and Conway calculated this standard deviation assum- 


ing that the distributions of Cc, and C. are Poisson in nature and 


therefore that the variances are expressible as: 
C 


Lee) le (8a,b 
and @(C,) = = ) 


Dee A 
TON rs 
where t is the duration of the measurement. Therefore they calcu-= 


lated  (c,) from the relation 





~ (C.) s EES pine ing Se (9) 


If the equipment used in the experiment were perfectly stable, with 
no changes in 7 due to changes in electronic characteristics, then 
all fluctuations in c, and a would be consistent with Poisson dis- 
tributions and these equations would apply. As pointed out in 
Section B, o°(C,) when calculated from equation 8a is too small. 

To calculate a realistic variance of C,; one must provide not only 
for the statistical variations due to the random nature of the radio- 
active decay process, but also for the variations in the resolving 
time of the coincidence unit due to changes in its electronic oper- 


ating characteristics. 
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If electronic instabilities were present and nuclear fluc- 
tuations were small in comparison thereto, then the variance of 


C, would be described adequately by the relation 


Cee Gale a 
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where Ceres denotes the standard deviation due solely to the elec- 
tronic instability of the coincidence unit. 

Since under experimental conditionsneither of the two limiting 
conditions obtains; i.e., perfect equipment stability in which case 
the variance of C, would be Poisson, or infinite counting time in 
which case nuclear fluctuations would be insignificant in comparison 
with OF ie (C,)s both a deviations must be taken into 
account in calculating ¢ (C,)- While the technique for the rigorous 
Statistical combination of these variances is not certain in spite 
of efforts to arrive at a valid analytical expression, the rather 
naive approach of adding them would seem to be a conservative ap- 


proach which applies to the experimental situation. Thus the ex- 





pression used in this report for calculating ia, ) is: 


@ é 2 
v(C,) = Semee CAae Telee 1) Ta (11) 
t t a 


It can be seen that this expression gives the correct standard 
deviation in the two limiting cases discussed earlier, and it is 
assumed that the error resulting from applying this expression in 


thevactual situation is "smai. 
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On the basis of the experimental determination of 7 , the third 
term of this expression is of the same order as the first two; hence, 
failing to take it into account results in unrealistically small esti- 


mates of the uncertainties in the experiment. 


Preparation and Check of New Computer Programs 
Conway prepared two programs in 1960-61 in 1604 computer assem- 
bly routine language. The purpose of one was to calculate the solid 
angle correction following the scheme of Rose [3 | to apply to the 
theoretical correlation function to account for the attenuation re- 
sulting from the finite dimensions of the detectors. The second 
program calculated the coefficients for a least squares fit [7 | of 
even order Legendre polynomials to W694) > 
As mentioned previously, difficulties were experienced in at- 
tempting to run these programs. In addition, the least squares pro- 
gram required the following prior manual computations: 
1. application of scaler dead time corrections 
2. correction of counting rates for background 


3. calculation of C, and T7,) 


&. calculation of C. and (Ca) 
2 
5. calculation of W, (0,) and Ho 95)] 
6. calculation of a weighted mean of W694) at each angle 


7. calculation of the standard deviation of the weighted 


mean at each angle 


The values of the weighted mean of W698 a? and the standard deviation 
of the weighted mean at each angle were the entering arguments for 
Conway's least squares program. 

The original intent was to program the above listed preliminary 
computations only, and to append Conway's least squares program as a 
subroutine. It was not possible to append his program to the FORTRAN 
program, however, due to the incompatibility of the two computer 


languages. A FORTRAN program was written to accept as input the total 
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counts registered by the scalers and thus save considerable time 
and achieve greater accuracy by avoiding the errors likely to creep 
into the lengthy manual computations. 

Separate programs were written which would compute the coeffici- 
ents for the even power cosine series as well as the coefficients for 
the even order Legendre polynomial series representation of W,, (9) > 
calculate the solid angle corrections to be applied to W469) » and 
calculate7 from the data obtained during coincidence resolving 
time runs. 

The above programs provided independent checks of Conway's 
least squares and solid angle correction programs. The new least 
squares program proved extremely useful during the reanalysis of 
earlier data as it permitted repeated reduction of the same experi- 
mental data with incrementally varied input parameters. 

The FORTRAN programs together with a glossary of terms is in- 
cluded in the Appendix(A thru E.) 

To check the operation of the FORTRAN programs, they were used 
to perform the same least squares and solid angle correction computa- 
tions as performed earlier by Verser and Conway, and the results com- 
pared with their earlier results, as follows. 

1. The same geometrical parameters as used by Verser and 
Conway were used as the entering arguments for the 
FORTRAN solid angle correction program. Good agree= 
ment was reached with Verser's hand-calculated values 
and Conway's computed values. 

2. The raw data (scaler counts) or Verser were used as 
the entering arguments for the FORTRAN least squares 
program. Good agreement was reached with Verser's 
values of the mean values of W 695) and with his 
coefficients for the Legendre series representation 
of W609) thus indicating proper operation of the 


preliminary and subsequent portions of this program. 
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3. The results of analysis of Conway's data with the 
FORTRAN program did not agree with the results of 
analysis which he reported. Some significant 
arithmetical errors were uncovered in the hand 
calculations which he was required to perform 
prior to entering his program. When these errors 
were corrected, agreement was realized in the re-= 
sults of analysis by the two programs, providing 
another check of the FORTRAN least squares program. 
Reanalysis of Earlier Data 
Once assured that the new FORTRAN programs were operating pro= 
perly, Verser's data and a portion of Conway's data were reanalyzed, 
including in the calculation of the variances the effect of the 
electronic instability of the coincidence unit. 
In checking over the data collected by Conway, it was noted that 
they were collected during two relatively short periods of time spaced 
nearly four months apart. During these four monthsa significant change 


in the strength of the Cou source (T ~5.2 years) occurred; this 


had not been corrected for in his ee ots. The reanalysis of 
his experiment was done using only his later set of data. 

It became apparent that the scatter of Conway's data about W,.(®) 
was greater than the standard deviations would imply. Consequently, 
estimates of the values of r{W(9;)| were calculated based upon the 
actual residuals between the experimental points and W,. 09) for his 


later set of data. This calculation, which took into account the elec= 


tronic variations of the equipment, resulted in values for 


o- | Me (,)| 
We 


approximately three times as large as the values which Conway quoted. 

The results of the reanalysis of the previous data are tabulated 
below and are presented graphically in Figure 6a and 6b. It is evident 
that the anomaly in the 12021350 region still exists in Verser's results, 
though it is not as significant statistically as previously reported due 
to the larger standard deviations. The analysis of Conway's data now 


shows no evidence of a dip. 
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Table 1 


Results of Reanalysis of Earlier Data 


Legendre Coefficients 
2 


Oo 4 


1.0000 + 0.0042 | 0.0960 + 0.0068 | 0.0337 + 0.0087 


0.162 
0.171 + 0.026 


The possibility was investigated that errors could have been intro- 













duced into Verser's data analysis through the use of incorrect parameters. 
In investigating this possibility, raw data gathered by both Verser and 

Conway were used as the entering arguments for the FORTRAN least squares 
Program, and these data were reduced while systematically varying the ex- 
perimental parameters to see what effect this would have on Wi, 69)- The 


results of this analysis are as follows: 


Table 2 
Effect on W,. (9) of Varying Input Parameters 


Parameter Varied % Variation in > Resulting 4% Variation 


Parameter 





It can be seen from the above that the shape of W,, 69) is relatively in- 


sensitive to changes in background and changes in the resolving time 
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of the scalers. However the shape of W,, 6%) varies appreciably 

when 7 varies. The reason for these changes in W,, 69) is clear 

from the following. If the numerical value used for 7 is greater 
than the actual 7 of the coincidence unit during the experiment, 
then too many accidental coincidences are subtracted from the total 
coincidences in calculating the number of genuine coincidences. 

Thus the apparent value of Co is smaller than was actually realized 


experimentally. Since the error introduced into C, is constant at 


A 
all angles, assuming that Ny and N, are independent of 9, which is 
approximately true, the error introduced into CC. is also constant 


G 
at all angles. Since in these experiments, the ratio of the acci- 


dental coincidence counting rate to the genuine coincidence counting 


rate is about 2:3 and it can be seen from the expression 


W, (180) 


Reh = W, (90) cia Reema Wy od (3) 


that W, (99) is more sensitive than W, 6189) to changes in the genuine 
coincidence counting rate, it is apparent that the calculated anisotropy 
increases sinometiie coincidence resolving time is increased. [t is 
stressed that the shape of W, (@) is distorted by using a value of ‘e 

in the data reduction which does not correspond to the value of 16 
actually characteristic of the coincidence unit during the experi- 

ment. The shape of W169) is independent of 7 as long as the true 
experimental value of 7 is used in calculating CA: 

The effect upon the shape of W,, 69) of using incorrect para- 
meters for the geometry of the experiment was investigated also, by 
progressively varying the entering arguments for the FORTRAN solid 
angle correction program. The results of this analysis are given in 


in Table 3. 
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Table 3 


Effect on Solid Angle Corrections of Varying Input Parameters 








Parameter Varied Amount Varied | Relative Change in 
| Absorption Coefficients 
| Q, Q, 
Linear absorption 
coefficient of the + 10% 0.00% 0.00% 
crystal 





Distance from source 0. 00% 0.00%, 
to crystal = 


Crystal radius + 1.0% 0.00% 0.00% 


Crystal thickness + 1.0% 0.00% 0.007 


It can be seen that the shape of W, 69) is relatively insensitive to 
changes in these geometrical parameters over the range investigated. 

It is concluded from these investigations of the effects upon the 
least squares analysis and solid angle correction computations of varying 
the entering arguments, that none of these can account for the dip re- 


ported by Verser. 
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H. Analysis of New Data 


The results of the analysis of the data obtained with the old 
slow coincidence unit in experiments in two separate quadrants are 
indicated graphically in Fig. 7. When the data from each of these 
two experiments were reduced separately, W, , (9) agreed with W169 
Within the statistical error of the experiment. Since the angular 
correlation function is symmetrical about 180°, it was possible to 
combine the results of the analysis of the data obtained during the 
experiment in the 90 - 180° quadrant with those obtained during the 
the experiment in the 180 - 270° quadrant. The purpose of combining 
the results in this way was to benefit from the decrease in r{w,.(0,)| 
resulting from the analysis of a greater quantity of data, thus being 


able to compare more precisely W, , (9) and We, 69) 


Due to the time interval between the two experiments, during 
which the source decayed significantly, the raw data from the two 
experiments cannot be combined. Consequently, the data obtained from 
each experiment were analyzed independently, then the normalized 
values of the weighted mean of W694? from the two experiments were 
combined using as weighting factors the reciprocals of their variances. 
These composite values of W569.) were then used as the entering 
arguments for the least squares subroutine to generate a new W, , (9) 


The results of the new experiments are tabulated below. 


Table 4 
Results of Analysis of New Data 


: Experimental Theoretical 


Quadrant 


We 1s0°-2707 


—_— =-———» A BE lemme e  E TF  e a HE 


| 0.151 + 0.029 0.165 +. OWS =) 0.158 4 008s 0-162 





Composite 






ITD ae ee a RN amet 


Anisotropy 








Legendre 
Coefficients 





1.0000 + 0.0057 | 1.0000 + 0.0040} 1.0000 
0.1062 + 0.0094 | 0.0987 + 0.0069; 0.0996 


-0.005 + 90.012 | 0.0036 0.0085] 0.0084 


I+ 
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Results of Analysis of New Data 


Figgas 


24 








Comparison of Experimental Results with Theory 


All of Verser's data, and the portion of Conway's data used in 
the reanalysis of his data (see Section G) were combined with the 
data obtained during the two new experiments according to the scheme 
discussed in the previous section. That is, the data from each of 
the four experiments were analyzed independently, then the weighted 
means of the normalized mean values of W(85) from the four experi- 
ments were used as the entering arguments for the least squares sub- 
routine. The results of this combination of these data are indicated 
graphically in fig. 8, and are tabulated in Table 5. No significant 
deviation from We 69) is evident. 

While the real demonstration of the proper functioning of the 
experimental equipment lies in the results of the new experiments in 
two quadrants, it is of interest that no significant deviation from 
W 1,69) occurs even when Verser's and Conway's results are combined 
with these. It is concluded that the dip which they reported is 
strictly a statistical fluke and that there is no actual physical 


effect involved. 


Table 5 


Results of Composite Data Analysis 
(Composite Results based upon 4 experiments) 


| | ee aes. 
| Experimental Theoretical 


nn LOA LL TEC Ns 
ee ee er a ty Se ee ee i EES ON ee NE VE oA cae ee EN a CTE A = Ry tees 9p ml Sem Ree ee vg ee eH - =. 
' 1 { 
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j 
Anisotropy 0.169 + 0.018 — o1162 | 
: Legendre | an 
| Coefficients | 
| a ) 1.0000 + 0.0028 1.0000 | 
| a, 0.1000 + 0.0047 0.0996 
a, 0.0183 + 0.0058 0.0084 


ee ge ne ee 
RE See 











Results of Composite Data Analysis 


igi on 
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IIil 
DEVELOPMENT OF NEW INSTRUMENTATION 
New Slow Coincidence Unit 


While investigating the earlier reported anomaly in W, (9) 7 
the 105-135" region, work proceeded on building and testing a new 
slow coincidence unit incorporating certain portions of the design 
reported by Garg | 131 . This design features a variable resolving 
time which is independent of the counting rate. The slow coinci-+- 
dence unit used by Verser and Conway lacked both of these features. 

The coincidence unit which was constructed differed markedly 
from Garg's only in that the separate channel amplitude discrimi- 
nators of Garg's unit were deleted since differential discrimina~ 
tion was performed earlier (Fig. 9 ). Figure 10 gives the com- 
plete circuit diagram for the unit which was constructed. 

The new slow unit was tested by incorporating it into the 
circuitry used for the experiment in the 90-180° quadrant in 
parallel with the slow coincidence unit which Verser and Conway had 
used. The correlation function measured by the two slow coinci= 
dence units agreed within the uncertainties in their resolving times. 
The results of numerous resolving time tests of the new slow coinci- 
dence unit indicated a greater long time drift than is characteristic 
of the old slow unit. This increased drift of 7 necessitates more 
frequent determinations of its magnitude and required using different 
values for T during the angular correlation runs, rather than using 
a mean value of 7 based upon a large number of resolving time runs 


such as was done with the old slow coincidence unit. 


Requirement for a Fast Coincidence Unit 
iv, General 
A coincidence unit of significantly shorter resolving time 
than that of the unit used by Verser and Conway was required to 
improve the statistical error in the experimental estimate of 
the correlation function, Wa 69.)- Better statistics were re= 


quired for the reinvestigation of the earlier reported anomaly 
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in the value of W695) for Nee in the 105+135° region, and 
for future experiments with short-lived nuclides produced in 
the US Naval Postgraduate School Van de Graaff generator. 

It was not possible to reduce the resolving time of the 
slow coincidence unit below a value ee which is determined 
by the width of the window of the differential pulse height 
analyzer which precedes it. This time corresponds to the dif- 
ference in time between the crossing of the lower voltage dis- 
crimination level of a pulse of the lowest permissable amplitude 
in one channel and the time of crossing of the lower discrimi- 
nator level of a simultaneous pulse of the highest permissable 
voltage in the second channel. Because of this irreducible 
minimum value for 7 » and since a resolving time larger than 

tbe: gives too many accidental coincidences for an acceptable 
variance of the experimentally-determined correlation function, 
a different mode of operation was required. The scheme used 
consisted of adding a fast coincidence unit in parallel with 
the slow unit and mixing the output pulses from the two units; 
this scheme reduced the accidental coincidence counting rate 
because of the decreased 7 of the fast unit, yet retained the 
amplitude: discriminating characteristics of the differential 
pulse height analyzer employed with the slow unit (Figure ll .) 

The fast coincidence unit was constructed after a design 
by Draper and Fleischer (14) and was adjusted for a resolving 
time of approximately 20 nanoseconds, thus decreasing cr by 
about a factor of 10 from that obtained with the slow units. 
Relation betweeng/W,,(0,)| and g(C,) 


It is desired to minimize 
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and N, are very small compared to the fractional deviation of 
Cus 
y eet ye CG (12) 
Wa Oey C 


Effect ono“ (6,) of Decreasing the Coincidence Resolving Time 


The variance of C is calculated from equation i 


G 
oC Ne Cr - aed! C Z " oA 
Get Fe A t 
It is evident from the relation C, = 27N.N iN, [equation 6 that 


decreasing the resolving time of the coincidence unit (T] by an 
order of magnitude results in a similar decrease in Cy: This 
decrease is reflected as a decrease in the third term of equation 
ll by an order of magnitude since t is constant in these experi- 


ments. ECS no change in the fractional resolving time fluctu= 
eled 


r 


by a factor of about 100. The genuine coincidence counting rate 


ation, e ), the second term of equation l11lwould decrease 
is independent of the coincidence unit resolving time, as long 
as the latter exceeds the average interval between the emissions 
of the two quanta; this is the case for all three coincidence 
units discussed in this report. Thus a decrease in C. also ree 


A 


sults in a decrease in C, since C, = oF mh ce This decrease in 


C,. results in a decrease in Cy, Therefore all three terms of 
Ct 

equation 11 decrease when 7 decreases. 

Relative Effect of Coincidence Resolving Time Drift 


The previous section indicated how (6) is decreased by 





decreasing T . Then from the relation 
o{W,,.(0,)| gC) 
(oi = G » equation 12, 
ex’ i G 
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it is evident, since C, is constant, that 


o{M<(0;)| 
ota? 


Ue le 


7 


did not vary was made in the last section to simplify the dis- 


has been decreased appreciably. The assumption that 


cussion. It is evident, however, that due to the large decreases 
in the values of all other terms in equation 11, including the 


6 
See , larger values of the 


ce term which multiplies 
latter can be tolerated while still decreasing very significantly 

the fractional variance of We: Thus the electronic stability 
of the fast coincidence unit is not as critical as it would be 


with a unit of larger coincidence resolving time. 


Operation of the Fast Coincidence Unit 

The basic theory of operation of the fast unit and the circuitry 
for the Draper:-Fleischer unit are given in detail elsewhere [14] and 
will not be repeated here. Either single channel becomes operative 
when the input pulse reaches 68% of its maximum amplitude, regardless 
of the value of that amplitude within the range 20-100 volts. The cire 
cuit is so designed that resolving times on the order of 100 times 
shorter than the rise time of the input pulses can be achieved. This 
characteristic is required since the rise time of the pulses from the 
sodium iodide scintillation crystals is on the order of 200-300 nano- 
seconds. 

Only the coincidences of 1.17 MEV and 1.33 MEV gamma»rays are of 
interest in this experiment with ee) Amplitude discrimination which 
was performed by the differential pulse height analyzers which preceded 
the slow coincidence unit selected the pulses corresponding to the 
photopeaks of the two gamma rays. All coincidences, regardless of the 
pulse height, were detected by the fast unit since no amplitude dis- 
crimination preceded it. The output pulses of the slow coincidence 
unit and those of the fast unit were passed through a mixer circuit 


which effectively functioned as another coincidence unit. An output 


* A block diagram of the fast unit is provided in Fig. 12. 
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pulse was obtained from the mixer only if there was, in a time interval 
of about 35 microseconds, a 20nanosecond coincidence between pulses 

of any height as ascertained by the fast unit and a 200-nanosecond 
coincidence between pulses corresponding to the photopeaks of interest 
as determined by the slow unit. Thus by combining these two units 

in this manner, an energy~selective counting system with a 20-nano- 
second coincidence resolving time which was independent of the counting 
rate was obtained. 

Certain modifications were necessary to the fast unit and are 
given in Figure 13. The last two stages of the Draper-Fleischer 
unit could not be made to function properly. This portion of the 
circuitry was therefore taken out, and a cathode follower was sub= 
stituted. 

The mixer circuit which contained circuitry for amplitude dis= 
criminating, shaping, amplifying, and mixing the fast coincidence 
pulse with the slow coincidence pulse was built up and tested stage- 
by-stage. The discriminator of the mixer unit is a uni-stable multi- 
vibrator with Zener diode stabilization. The inverter stage in the 
fast coincidence channel of the mixer results in a negative pulse 
of approximately 35-microsecond duration, the duration varying with 
the discriminator setting. This pulse is then mixed with a 0.7 microe- 
second negative pulse from the slow coincidence unit. Discrimination 
of the mixed output is performed by the scaler counting these fast- 


slow coincidences. The design of the mixer circuit is given in Fig. 14. 


Testing and Alignment of the Fast Coincidence Unit 

The alignment of the fast coincidence unit consisted of two 
major portions-achieving linearity in each channel, and assuring the 
detection of true coincidences by the coincidence stage. 

Draper and Fleischer | 14 describe a means of testing for 
linearity; i.e., ensuring that the fast discriminator fires when the 
input reaches the same predetermined fraction of its full amplitude, 
over a wide range of amplitudes. They prescribe a mercury-relay pulser 
delivering a five-volt negative step with a rise time of seven nano-= 
seconds which is split, one portion triggering a 10-mc ringing cir- 


cuit which provides the time reference, the other going through an 
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attenuator, linear pulse amplifier, and the coincidence unit 
(Figure 15 -a)., 

Several modifications of this test scheme were tried, due to 
the nonavailability of a pulse generator which could deliver a pulse 
of sufficiently large amplitude and short rise time to trigger the 
ringing circuit satisfactorily, and due to the lack of a compensated 
attenuator. These modifications met with varying degrees of success. 
When a double pulse generator was used as a signal source, and the 
amplitudes of the inputs to the separate channels of the fast unit 
were varied by adjusting the gains of the linear amplifiers, dis- 
tortion of the wave shapes resulted when the coarse gain settings 
were varied. This made it impractical to vary the amplitudes of 
the coincidence unit inputs by varying the gains of the linear 
amplifiers since the time reference for determining the delay within 
the coincidence unit itself was thus lost. Still greater wave shape 
distortion resulted from varying the gain of the double pulser. 

The procedure which has finally emerged as workable is outlined 
in Figure 1l5-b and consists of triggering the oscilloscope with 
the output of the fast unit discriminators and viewing the relative 
delay of nuclear pulses of varied amplitudes. The problem of wave 
shape distortion when the coarse gain settings of the linear ampli- 
fiers were adjusted has been avoided by limiting the adjustment of 
delay-equivalence to the 10:volt range of amplitudes used experi- 
mentally. Since a 10-volt window is used in the differential pulse 
height analyzer experimentally, an equivalent time delay over this 
range is all that is required. The inputs to the fast unit are varied 
over this 10-volt range by adjusting only the fine gain controls, 
avoiding the distortion of wave shape which accompanies changes of the 
coarse gain settings. The adjustment of the fine gain controls can 
be avoided by providing a nuclear source or a combination of sources 
with different decay energies which results in different amplitude 
signals being received by the fast unit simultaneously. The relative 
delay introduced into these pulses of different amplitudes is then 


observable directly. 
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An alternate means of accomplishing this alignment consists of 
using test pulses from a double pulser, triggering the oscilloscope 
with the trigger output pulse of the double pulser and viewing the 
output of the fast discriminator of the fast coincidence unit in 
each channel. Again, the input amplitude is varied by changing the 
settings of the fine gain controls on the linear amplifiers. 

The procedures already described ensure only that the delay 
which a pulse experiences as it passes through the time compensating 
circuit of one channel is independent of amplitude. The other major 
alignment requirement, which cannot be met until each channel has 
been aligned separately, is the assurance that the true coincidences 
are being detected by the unit. This requires identical total delay 
from the time of emission from the source through all the circuitry 
associated with each of the two channels and through the fast coinci- 
dence unit itself up to the coincidence stage. 

To determine the ability of the fast unit to detect true coinci-= 
dences before it was completely aligned, it was used in parallel 
with a slow unit in the angular correlation experiment in the 90-180° 
quadrant. While it generally tracked the slow unit, it did so errati- 
cally. 

Variable delay lines have been installed to facilitate the align- 
ment for equivalent delay in the two signal channels to the coincidence 
stage of the fast unit. The alignment is accomplished using the experi- 
mental coo source with the experimental geometry and normalizing the 
total coincidence counting rate measured by the fast-slow coincidence 
detecting scheme to the genuine coincidence counting rate determined 
by the slow coincidence unit. After normalizing the total coincidence 
counting rate of the fast~slow combination in this manner, and making 
a graph of this versus the setting of the variable delay line in one 
of the time compensation channels of the fast unit, a plateau is ob- 
tained the width of which corresponds to the resolving time of the 
fast coincidence unit. The delay line is then adjusted to the center 


of this plateau as determined from the graph. 
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Thus procedures have been developed for aligning the time com- 
pensation circuit of each channel for linearity by means of either 
a nuclear source or a double pulser, and for aligning the two chan-= 
nels of the fast unit for equivalent delay by means of a nuclear 


Source. 


Equipment Status 

The experimental equipment described herein including the new 
fast~slow coincidence detecting scheme is functioning well enough 
to be used in the experiments for which it was assembled; namely, 
the angular correlation of short: lived nuclides activated in the 
Van de Graaff generator. Due to limited unavoidable drifts in the 
electronic equipment, frequent checks must be made of operating 
characteristics, and appropriate adjustments made as necessary to 
ensure the closest possible control over the experimental environ= 
ment. The checkout procedures which have proved extremely useful 
in ensuring that easily~overlooked actions are not overlooked and 
that checks are made in the same sequence for the greatest possible 
reproducibility of procedures, should be continued. 


Figure l6illustrates the layout of the major electronic portion 


of the experimental equipment, all rack-mounted and mobile to facili- 


tate its use with or without the Van de Graaff generator. 
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16 Electronic System 


Fig. 


FORTRAN PROGRAM 


Appendix Part A 


Least Squares Analysis 
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Results of Analysis by the Least Squares Program 
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1. Check of input p 


TAUC . 
STAUC 
TAUD1L 
TAUD2 
BACK] 
BACK2 
T 

2. Check data input 


7 


MD 
CT 


3. Results of analy 
NDOT 1 PRIME 


NDOT 2PRIME 
CTOTALDOT 
CACCDOT 


CGENDOT 


S IGMAW 


WBAR 


Appendix E FORTRAN Program 
Glossary of Terms for the 
Least Squares Analysis 


arameters 


oincidence unit resolving time (sec) 
Standard deviation of 7 (sec) 

System dead time, Channel 1 (sec) 
System dead time, Channel 2 (sec) 
Background, channel 1 (counts per sec) 
Background, channel 2 (counts per sec) 
Duration of run (sec) 


of N,N,and CT and input parameters 


Total counts, channel 1 
Total counts, channel 2 
Total coincidence counts 
sis 
Count rate, channel 1, corrected for system dead 


time and background (counts per sec) 
Corresponds to N, in text 


Same, channel 2 
Co 
oh 
CG 
G 
ss soe een 
di 
AS) 
G 
r{W, (0:9 = hed x [Rex] 
/ ee 
(0.))- 
on (6;)= as M 


a) 


| 
SIGMAWBAR c (0, | : : 
ex L 








N 
affw,c@p) we.) 


WBARB T |W, @,)] = N (N - 1) 








where N is the number of runs at the angle Q 
EXPERIMENTAL FUNCTION W, (9) normalized to W, , (99) and 
evaluated at 5° intervals. 


THEORETICAL FUNCTION W169) normalized to Wb?) and 


One 
evaluated at 5 intervals. 


NORMALIZED EXPERIMENTAL WBARS AND SIGMAS 


WBAR Wi (9? normalized to W699) 


ANISOTROPY - P cue 


=———— - 1 | based upon Legendre series 
'W, (90) 
= vis representation of Wis: 
ANISOTROPY - C 


of W 


Same, based on Cosine series representation 
rs. 


SIGMA Standard deviation of ANISOTROPY - P 
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